Coulomb Suppression of NMR Coherence Peak in 
Fullerene Superconductors 
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The suppressed NMR coherence peak in the fullerene superconductors is explained in terms of 
the dampings in the superconducting state induced by the Coulomb interaction between conduc- 
tion electrons. The Coulomb interaction, modelled in terms of the onsite Hubbard repulsion, is 
incorporated into the Eliashberg theory of superconductivity with its frequency dependence consid- 
ered self-consistently at all temperatures. The vertex correction is also included via the method of 
Nambu. The frequency dependent Coulomb interaction induces the substantial dampings in the su- 
perconducting state and, consequently, suppresses the anticipated NMR coherence peak of fullerene 
superconductors as found experimentally. 



It is generally accepted that the superconducting prop- 
erties of the fullerene superconductors can be understood 
in terms of phonon-mediated s-wave pairing The 
l/(TiT) for s-wave superconductors, where Tj is the nu- 
clear spin-lattice relaxation time and T is the tempera- 
ture, is expected to show a peak as T is lowered below 
the transition temperature, T c , and is constant above T c . 
It is referred to as coherence peak or Hebel-Slichter peak. 
The expected coherence peak is found substantially sup- 
pressed for fullerene superconductors (|-|^] . The underly- 
ing mechanism of the suppression, however, is not clearly 
understood yet. The present Letter addresses the prob- 
lem of NMR coherence peak suppression in fullerene su- 
perconductors and explains it in terms of the dampings 
in the superconducting state induced by the frequency 
dependent screened Coulomb interaction between con- 
duction electrons. 

The maximum of the normalized relaxation rate, 
R s /R n , in the limit of zero applied magnetic field was 
estimated to be 1.1 — 1.2 for fullerene superconductors 
HQ, where R = T-f 1 is the relaxation rate, and the sub- 
scripts s and n, respectively, refer to the superconducting 
and normal states. A similar behavior was also observed 
in /iSR (muon spin relaxation) experiments p| . An anal- 
ysis based on the phonon-mediated Eliashberg theory ||] 
gives T c /wph ~ 0.2, and the characteristic phonon fre- 
quency u>ph ~ 100 cmT 1 0. This seems inconsistent with 
the commonly held view of aj p h ~ 1000 cmT 1 for fullerene 
superconductors [|l|-Q, which implies {R s /R n )max ~ 2 — 3 
PL The increase of (R s / R n )max relative to its normal 
state value, 0.1 — 0.2, is therefore suppressed by an order 
of magnitude from the expected value. Closely related is 
the observation of the unexpected sensitivity of R s /R n to 
the applied magnetic field. The coherence peak is found 
completely suppressed at only around 5 T for Rb2CsCeo, 
which is at least an order of magnitude smaller than the 
expected value Q. It may be understood on the ground 
that the zero field coherence peak is substantially sup- 
pressed as mentioned above. A simple weighted average 
of the suppressed intrinsic rate from the superconducting 
region and the normal rate from the normal region inside 



the vortex core will give the experimentally observed field 
dependence of R s /R n [Q. 

In order to understand the coherence peak suppres- 
sion at zero magnetic field, let us list the possible fac- 
tors that are known to affect the NMR coherence peak 
of R s /R n pC|p| , and check if any of those can explain 
the suppression. The NMR coherence peak suppression 
may be attributed to (a) momentum anisotropy of super- 
conducting gap including non-s-wave pairing, (b) time 
reversal symmetry breaking such as magnetic impuri- 
ties or applied magnetic field, and/or (c) the damping 
effects in superconducting state. For fullerene supercon- 
ductors, the superconductivity is of the phonon-mediated 
s-wave pairing, and due to the orientational disorder of 
Ceo molecules, the fermi surface anisotropy is not strong 
enough to suppress the coherence peak JTTJ . The time 
reversal symmetry breaking can not explain the suppres- 
sion either because there are no magnetic impurities in 
the fullerene superconductors, and we are considering 
the case of zero applied magnetic field. Because either 
the gap anisotropy or time reversal symmetry breaking 
can not explain the coherence peak suppression in the 
fullerene superconductors, the near absence of the coher- 
ence peak should be due to the damping effects. 

The dampings of an electron come from the scatter- 
ings of the electron with the phonons, impurities, and/or 
other electrons. The damping from the scatterings of 
electrons with the phonons, that is, the electron- phonon 
interactions, is not strong enough: The dimensionless 
electron-phonon coupling constant, A, of the fullerenes 
is estimated to be 0.5 — 1 j| [|], far smaller than ~ 2 
needed to suppress the NMR coherence peak |L^ ]. The 
far infrared reflectivity measurements of DeGiorge et al. 
[p[ show that the ratio 2A /k B T c « 3.44 - 3.45 both 
for K 3 C6o and RbsCeo- It is very close to the BCS 
value of 3.52 and implies that A can not be as large as 
2. The impurity scatterings smear out the fermi surface 
anisotropy and can not suppress the coherence peak as 
explained previously. There remains, therefore, only one 
possibility for inducing the dampings required to sup- 
press the NMR coherence peak in the fullerene supcrcon- 
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ductors: The scatterings between electrons due to the 
Coulomb interactions. This idea is indeed verified in 
our detailed Eliashberg-Nambu (EN) calculations as will 
be detailed below. 

The NMR relaxation rate for a superconductor with a 
finite bandwidth of B is given by 



1 r , 9f F (e) J [ e8(e) 
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) is the fermi distribution func- 
e) = tan" 1 ( J B/2Z(e) A /e 2 - A(e) 2 ) , 



where /f(e) = 1/(1 ' 
tion, /? = l/fc B T, 0( 

and Z(u>) and A(w) are, respectively, the renormaliza- 
tion and gap functions. The finite conduction band- 
width with a constant density of states (DOS) is ex- 
plicitly considered through the factor of 0, which is tt/2 
for the usual case of infinite bandwidth superconductors. 
For fullerene superconductors, the fermi energy ep — 
8/2 k 0.2 — 0.3 eV and the average phonon frequency 
w p h w 0.05 — 0.15 eV. Consequently, 0Jph/£F ~ 1 for 
fullerenes unlike conventional metals where ujph/eF <C 1. 
When ujph/ep ~ 1, the phonon vertex correction becomes 
important because the Migdal theorem does not hold, 
and the Coulomb interaction should be considered more 
carefully because the validity of the Coulomb pseudo- 
potential, u* = fi/[l + (xki(ep/u> p h)], is unclear. In the 
present work concerned with the narrow bandwidth su- 
perconductor of fullerenes, therefore, the vertex correc- 
tion is incorporated into the Eliashberg theory via the 
method of Nambu |Q . The Coulomb interaction is mod- 
elled in terms of the onsite Hubbard repulsion and is in- 
cluded in the theory with its frequency dependence con- 
sidered self-consistently. As far as we are aware of, this 
is the first self-consistent calculation of 1/ (XiT) with the 
Eliashberg formalism including the Coulomb interaction 
and the vertex correction. The frequency dependence 
of the screened Coulomb interaction, which comes from 
the polarization diagrams of both the normal and pairing 
processes of renormalized electrons, is important because 
the electron-phonon and Coulomb interactions vary on a 
comparable frequency scale for a narrow bandwidth su- 
perconductor of UJph/^F ~ 1. 

The NMR coherence peak below T c is due to the in- 
creased DOS in the superconducting state, that is, due 
to the smallness of the denominator of Eq. ([j]) when 
e w A(e). In order for the peak to be suppressed, 
therefore, the vanishingly small denominator should be 
avoided. The \J e 2 — A(e) 2 of Eq. (|l|) may not vanish 
when there is a damping, that is, non-zero A2(e) for 
e w Ai(e), where Ai and A2 are, respectively, the real 
and imaginary parts of the gap function. The dampings 
in the superconducting state responsible for the NMR co- 
herence peak suppression is greatly increased when the 



frequency dependence of the screened Coulomb interac- 
tion is retained for ujph/tF ~ 1- There have been several 
works which emphasize the importance of including both 
the electron-phonon and electron-electron interactions in 
understanding the fullerenes [T^ ]. 

The EN equation can be written in the Matsubara fre- 
quency as: 
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) is the Matsubara frequency, 

On = tan" 1 (B/2Z n ^/p 2 n + A 2 ^j, and X ph (n - m) = 

f °° dfla 2 F(fl)2fl/ [n 2 + (p n - p m ) 2 ] is the pairing ker- 
nel due to the electron-phonon interaction. Eq. (|2|) is of 
the same form as the theory used to study the spin fluc- 
tuation effects on superconductivity Jl6| |. The X c h(n — m) 
and X sp (n — to) are, respectively, the interactions in the 
charge and spin channels due to the Hubbard repulsion, 
and are determined self-consistently as 

X ch (k) = UN F {1/2 -( Xn+Xs ) 

+ (Xn+Xs) 2 \n[l + l/(Xn+Xs)}} , 

X sp (k) = UN F {1/2 + ( Xn - Xs ) 

+ (X«-X,) 2 ln[l-1/(X«-X,)]}, (3) 

where Xn {k) and Xs (k) are the dimcnsionless susceptibil- 
ities from, respectively, the normal and pairing processes 
given by 
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The Z m on the right hand side of Eq. (0) represents the 
vertex correction of Nambu, which we take as its form in 
the normal state for simplicity ju| . Including the vertex 
correction enhances the transition temperature in accord 
with the previous works Jl7j . If we neglect the vertex cor- 
rection, the Z m should be put equal to 1. Self-consistent 
solution of Eq. (||) together with Eqs. (||) and (Q) gives 
Z(ip n ) and A(ip n ) in the imaginary frequency. To obtain 
Z(u>) and A(lo) in the real frequency, we perform the an- 
alytic continuations using the iterative method of mixed 
representations (l^Jl^]. It is more efficient than solving 
the Eliashberg equation directly in the real frequency. 
The details of the Eliashberg-Nambu formulation in the 
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imaginary frequency and its analytic continuation will be 
reported separately. 

The EN equation of Eq. ^ is solved self-consistently 
via iterations with a set of the phonon spectral func- 
tion a 2 F(£l), Hubbard repulsion U, and impurity scat- 
tering rate t _1 . To model the fullerene superconduc- 
tors, we take a 2 F(Q) = Y,l=i°Z F v(tyi where F »( n ) 
is the truncated Lorentzian centered at w v with the 
broadening T = uj,s/5, the cutoff frequency T c = 3T, 
and dSlF v (£l) — 1 plfl . Various theoretical and ex- 
perimental estimates do not agree well with each other 
in terms of distribution of coupling strength o? v among 
the different modes. These estimates show, however, 
that the phonon spectra derived from the intramolec- 
ular A g and H g modes are distributed over 0.03 — 0.2 



eV with the total A in the range 0.5 — 1 j2j |4|. In 
view of this, we represent the phonon modes with three 
groups centered around lo v = 0.04, 0.09, 0.18 eV, and 
2N F a 2 /uj u = 0.3 A s ,0.2 A s ,0.5 A s , respectively, for v — 
1, 2, 3. Note that Y^ v=1 2Npal/u) v = X s . This choice 
of a 2 F(fl) gives the logarithmically averaged phonon fre- 
quency u>i n w 0.094 eV, which is a representative value 
of the various estimates of u>i n - The X s is set to give 
T c ss 20 — 40 K for a given U. We take the fermi energy 
€p = B/2 = 0.25 eV in the present calculations, which 
gives the ratio u>i n /ep ~ 0.38. 

We put r _1 = for simplicity because the results are 
insensitive to the impurity scatterings. Even though the 
Anderson theorem does not hold exactly because of the fi- 
nite bandwidth and the Hubbard repulsion in the present 
theory, the thermodynamic properties are still insensitive 
to the impurity scatterings. We take UNp = 0.31 and 
A s = 0.71 in the numerical calculations reported below. 
We find T c = 0.0031 eV and 2A /k B T c = 3.1. It is inter- 
esting to note that the present theory gives rather small 
2Ao/fcsT c value which lies at the low end of the vari- 
ous estimates of the gap values. In the present study, 
the long wavelength contribution of the Coulomb inter- 
action is not considered explicitly. Therefore, U should 
be taken as a screened value. The previous estimates give 
UNp w 0.3 -0.4 §. 

Fig. 1 shows the superconducting gap function A(uj) 
as a function of co at T = 0.001 eV obtained by solving 
the Eliashberg-Nambu equation. We took 220 Matsub- 
ara frequencies to solve Eq. (Q) by iterations and dis- 
regarded 20 high frequency data to avoid boundary ef- 
fects. The analytic continuations were carried out with 
2000 frequencies in the range between — 0.6 eV. Fig. 
1(a) is for an interacting system of UNp = 0.31 and 
A s = 0.71, and 1(b) is for a non-interacting case of U 
= and X s = 0.35, which gives T c = 0.0032 eV and 
2Ao/fcsT c = 3.7. The solid and dashed lines, respec- 
tively, represent the real and imaginary parts of A(u>). 
The three peaks in A(u) reflect the three peaks in the 
phonon spectral function a 2 F(fl). Note the difference 
in A2(w) between the two cases for e ss Ai(e): For the 
interacting case (a), the superconducting state has quite 



a strong damping of [A2/A1I S3 0.05 around e « Ai(e) 
due to the Coulomb interaction even at the low temper- 
ature of T/T c w 0.3, while A2/A1 = around e w A x (e) 
for the non-interacting case of (b) because the thermal 
fluctuations are quenched at this temperature. As T is 
increased, IA2/A1I is further increased due to the ther- 
mal fluctuations. Substantial A2 is what suppresses the 
NMR coherence peak as shown in Fig. 2. 
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FIG. 1. The superconducting gap function, A(uS), at 
T = 0.001 eV. The solid and dished lines, respectively, stand 
for the real and imagninary parts. We take UNf = 0.31 and 
A a = 0.71 for Fig. 1(a), and U = and A a = 0.35 for 1(b), so 
that the two cases have similar T c . The substantial damping 
of IA2/A1I ~ 0.05 around e ~ Ai(e) for the interacting case 
of (a) can be seen more clearly in the inset of Fig. 1(a). 

The A(lu) and Z(u>) obtained above are then used 
to calculate the relaxation rate T\ using Eq. ([!]). Wc 
show in Fig. 2 the normalized relaxation rate R s /R n as 
a function of the reduced temperature T/T c . The solid 
line is for UNp = 0.31 and X s — 0.71 corresponding 
to Fig. 1(a). As expected, the substantial damping in 
the superconducting state in the interacting system sup- 
presses the NMR coherence peak so that the maximum 
of R s /R n S3 1.15. By comparison, the corresponding 
curve for the non-interacting case of Fig. 1(b) (Z7 = 
and A s = 0.35) exhibits a much more pronounced co- 
herence peak as shown by the dotted curve, a hallmark 
of the weak-coupling s-wave BCS superconductors. It is 
clear that the strong Coulomb interaction can suppress 
the NMR coherence peak for phonon-mediated s-wave 
superconductors with a modest electron-phonon coupling 
constant. As U is increased, the NMR coherence peak is 
further suppressed. The dashed curve shows R s /R n for 
UN F = 0.4 and A s = 0.8, which gives T c sa 0.003 eV. We 
now show some of the available experimental data against 
the theoretical curves in Fig. 2. The filled circles and 
squares represent, respectively, the data by Stenger et al. 
from 13 C NMR of Rb 2 CsC 60 at B = 1.5 T and B = 3 T, 
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and the open up and down triangles represent the data 
by Sasaki et al. from 13 C NMR of K 3 C 60 at B = 2.93 T. 
It seems that the experimental data can be well described 
by U Np w 0.3 - 0.4 and A s w 0.7 - 0.8. We note that 
the Sasaki's data show somewhat more suppressed NMR 
coherence peak compared with the Stenger's data, while 
the Stenger's data show a bit narrow width compared 
with the theoretical calculations. 
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FIG. 2. The normalized relaxation rate, R a /R n , as a func- 
tion of the reduced temperature T/T c . The solid curve is 
for UNf = 0.31 and A s = 0.71 corresponding to Fig. 1(a), 
the dotted curve for non-interacting case of UNf = and 
A,, = 0.35 corresponding to Fig. 1(b), and the dashed curve 
is for UNf = 0.4 and X a — 0.8. The substantial dampings in 
the superconducting state for the interacting cases suppress 
the NMR coherence peak as can clearly be seen by comparing 
above curves. The filled circles and squares are the data from 
Stenger et al., and the open up and down triangles are those 
from Sasaki et al. See the text for more detailed discussions. 

We point out that Mazin et al. have analyzed the 
NMR data using the Eliashberg theory with two peaks in 
a 2 F(Q), one high frequency peak at around 1000 cm -1 
of corresponding A = 0.5 and another very low fre- 
quency peak around 40 cmT 1 of A = 2.7 ]19] . Such 
a low frequency mode with the strong coupling gives 
Wl n = 66 cmT 1 in sharp contrast with the commonly 
held view of u>i n ~ 1000 cm" 1 . The present work explains 
the NMR coherence peak suppression without such a low 
frequency mode, and confirms that the BCS-Eliashberg 
framework of uj p h ~ 1000 cm~ 1 and A ~ 0.5 — 1 can pro- 
vide a consistent description of the wide range of exper- 
imental data including T%, provided that the frequency 
dependent Coulomb interaction and the vertex correction 
are properly incorporated. 

To summarize: We have extended the standard Eliash- 
berg theory for narrow bandwidth superconductors by 
including the frequency dependent screened Coulomb in- 
teraction together with the electron-phonon interaction, 
and by including the vertex correction via the Nambu's 
method. We then solved the Eliashberg-Nambu equa- 



tion self-consistently at all temperatures to obtain the 
gap and renormalization functions, A(w) and Z(uj), re- 
spectively, which are used to calculate the nuclear spin- 
lattice relaxation rate, . The frequency dependent 
Coulomb interactions between conduction electrons in- 
duce the substantial dampings in the superconducting 
state and, consequently, suppress the NMR coherence 
peak in the fullerene superconductors. The present work, 
therefore, has shown that the T\ experiments can be 
understood with the view of uj p h ~ 1000 cmT 1 and 
A ~ 0.5 — 1. It remains to be seen if other experimental 
data can also be understood with the view. 
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